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1 . 0  INTRODUCTION 


The  radar  data  generator  calculates  range  and  azimuth  given 
the  geographic  coordinates  of  a  radar,  the  simulated  stereographic 
coordinates  (X,Y)  and  height  of  an  aircraft, and  the  geographic 
coordinates  of  the  JSS  region  center.  The  generation  of  a  radar 
report  involves  a  coordinate  transformation,  a  rotation  and  a 
conversion.  The  transformation  is  from  stereographic  coordinates 
into  a  rectangular  coordinate  system  whose  axes  are  directed  north, 
east  and  vertical  from  the  region  center.  The  rotation  maps  the 
coordinates  from  this  region  centered  rectangular  coordinate  system 
to  a  system  whose  axes  are  directed  north,  east  and  vertical  at 
the  radar  site.  The  conversion  from  Cartesian  to  polar  coordinates 
provides  range  and  azimuth.  This  report  derives  the  transformation 
and  rotation  formulary  to  be  used  for  radar  data  generation  in  the 
JSS. 

Background 

The  algorithm  used  for  generr,  of  simulated  radar  data  is 
in  effect  the  inve  se  of  the  coordinate  conversion  and  transforma- 
_  used  to  project  radar  reports  onto  a  common  coordinate  plane. 

Radar  reports  from  a  network  of  sensors  are  typically  converted 
to  common  plane  coordinates  by  stereographic  projection.  The  equa¬ 
tions  for  determining  locations  on  a  stereographic  plane  were 
developed  for  a  sphere  because  the  mathematics  were  simpler  than 
for  an  ellipsoidal  earth  model.  The  equations  for  stereographically 
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projecting  a  point  referenced  to  the  earth's  surface  reflect  a 
two-step  process:  (a)  points  on  the  ellipsoid  are  conformally 
mapped  onto  a  sphere,  called  the  conformal  sphere,  and  then  (b)  the 
mapped  points  on  the  conformal  sphere  are  stereographlcally  pro¬ 
jected  onto  a  plane  tangent  to  the  conformal  sphere.  The  radius 
of  the  conformal  sphere  is  calculated  to  minimize  distortion  in 
the  projection  process.  This  radius  is  shown  in  another  document  to 
very  inversely  with  region  size. 


2.0  COORDINATE  TRANSFORMATION 


Assume  that  one  is  given  the  position  (X,Y)  of  an  object 
stereographically  projected  on  a  plane  tangent  to  a  conformal 
sphere  of  radius  Eo>  The  problem  is  to  find  the  range  and  azimuth 
coordinates  (P,0)  of  the  object  measured  with  respect  to  a  radar 
on  the  surface  of  the  earth. 

Figure  1  depicts  the  stereographic  projection  geometry.  The 
ROCC  region  center  is  the  geographic  location  of  the  center  of  the 
smallest  circle  which  encompasses  all  radar  sensors  tied  into  the 
ROCC.  The  center  of  the  tangent  plane  is  the  ROCC  region  center, 
point  A.  The  planar  location  of  an  aircraft  at  point  T  in  Figure  1 
is  described  by  the  vector  AP  and  is  determined  as  a  function  of 
the  angle  i|<  and  conformal  sphere  radius  as  follows: 

AP  =  2E  Tan  (4*/2)  (1) 

o 

where  point  P  is  the  aircraft's  planar  location.  To  describe  point 
P  in  vector  notation  with  respect  to  a  rectangular  coordinate  sys¬ 
tem  centered  at  0,  the  ENZ  coordinate  system  depicted  in  Figure  2 
is  introduced.  Let  AP  be  defined  in  terms  of  common  plane  coor¬ 
dinates  (X,Y)  as  follows: 

AP-XE+YN+OZ 


AP 


X 

Y 

0 
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. . 


where  Z  is  normal  to  the  surface  of  the  sphere  at  the  region  center 
and  parallel  to  OA,  and  the  EN  plane  is  parallel  to  the  tangent  plane 
with  the  N  axis  directed  due  north. 

Figure  2  shows  the  vector  OP  in  the  rectangular  coordinate  system. 
OP  is  found  by  passing  a  line  through  point  0  and  point  B  and  inter¬ 
secting  with  the  tangent  plane  at  point  P.  The  vector  OP  may  be 
described  by: 

OP=xF+YN  +  |  AO  |  Z 

where: 


I  AO  I  =  2E 
'  1  o 


To  simplify  subsequent  analysis,  all  coordinates  will  be  norma¬ 
lized  and  referenced  to  a  conformal  sphere  of  unity  radius.  This  is 
accomplished  by  dividing  the  scalar  components  of  the  vectors  by  the 

radius  of  the  conformal  sphere  E  .  Thus  the  normalized  vectors  OP., 

o  N 

and  AP^  are  defined  by: 


U  -  X/E 

o 

V  =  Y/E 

o 

W  =  (2E)/En 
o  o 


2 


(2) 
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and : 


The  vector  CT  in  Figure  1  is  constructed  from  the  center  of 
the  conformal  sphere  to  an  aircraft  at  point  T.  The  magnitude  of 
CT  is  the  sum  of  the  aircraft  height  (H^)  and  the  radius  of  the 
earth  as  shown  in  Figure  3.  Normalizing  the  vector  CT  as  described 
previously,  the  magnitude  of  | CT^ ]  is  as  follows: 

|CTnI  =  (Eg  +  y/Eo  (4) 


where : 

E  =  Earth  radius, 
s 

H^,  *  Aircraft  height,  and 

E  •*  Conformal  sphere  radius 
o 

The  vector  CT„  may  be  defined  as: 

N 


IcthIcbn 


where  CB  is  a  unit  vector  parallel  to  CT  and  where: 


(5) 
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Let  us  now  derive  the  scalar  components  ,  C^,  and  of  the 
vector  CB^  in  terms  of  CO^  and  OB^  as  shown  in  Figure  4: 

Vectors  OB^  and  0PN  are  in  the  same  direction  but  with  differ¬ 
ent  magnitudes  so  OB.,  can  be  defined  as: 

N 


°Bn  -  kopn 


K 


u 

V 

2 


Thus: 


Ki 

lo>'Nl 


(b) 


To  solve  for  K  we  note  that  the  triangles  ABAO  and  AAPO  in  Figure  4 
are  similar.  The  triangles  share  angle  AOB,  and  angles  OAP  and 
OBA  are  both  right  angles;  thus: 


I  °°N '  .  inJ 

Ki  isn' 

From  (6)  and  (7) 

15J  15J2 

^  m  .  "  _  2 

I  I  I  I  *“ 


(7) 


13 


By  the  law  of  right  triangles: 


lorj2  -  Ioa/  +  |apn|2 


K'2 

l%l2  +  l“»l2 


From  Equation  3: 


I apn | 2  -  U2  +  v2  +  o2 


and  vector  OA^  is  the  diameter  of  a  unit  radius  circle  so: 


\0\\2  -  * 


Thus,  K  is  found  to  be: 


2  2 
4  +  U  +  V 


From  equation  (6)  OB^  is  defined  as: 


4U/(4  +  U2  +  V2) 
4V/(4  +  U,  +  V,) 
8/(4  +  U  +  V  ) 
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From  Figure  4,  we  note  that: 


where: 


CB»  ■  obn  *  % 


Thus, 


cbn 


4U/(4  +  ul  +  V2) 

4V/(4  +  IT  +  V% 

(4-U  -V  ;/ (4  +  U  +  V  ) 


(11) 


Coordinate  Rotation 

Having  obtained  the  unit  vector  CB^  that  is  directed  toward  the 
aircraft  in  the  normalized  ENZ  coordinate  system,  the  slant  range 
and  azimuth  (p»0)  from  a  radar  site  to  the  aircraft  can  be  calculated. 
The  ENZ  coordinate  system  is  centered  at  the  ROCC  Region  Center;  its 
axes  are  directed  toward  east,  north  and  normal  to  the  surface  of 
the  sphere  at  the  region  center.  The  latitude  and  longitude  of  the 
region  center  are  defined  as  (^o»0o)  respectively. 

To  calculate  (p,0),  the  ENZ  coordinate  system  must  be  rotated  into 
a  system  whose  axes  are  directed  east,  north  and  normal  to  the  surface 
of  the  sphere  at  the  radar  sites.  are  defined  as  the  latitude 

and  longitude  of  the  radar  site. 


lb 


Positive  longitudes,  X,  will  be  measured  counterclockwise,  or 
west,  from  the  prime  meridian,  X-0°  located  in  Greenwich,  England. 
Positive  latitudes,  0,  will  be  measured  north  from  the  equator,  0-0. 

Figure  5  depicts  a  set  of  rotations  that  align  the  ENZ  axes 
toward  east,  north  and  normal  to  the  surface  of  the  sphere  at  the 
radar  site. 

In  summary  the  sequence  of  rotations  is: 

1)  CCW  about  E,  0  ,  to  align  N,  with  the  earth's  North  axis 

o 

2)  CCW  about  N^,  XQ,  so  that  Z^  lies  on  the  prime  meridian 

3)  CW  about  N„,  \  ,  so  that  Z,  lies  at  longitude  \ 

2  s  i  s 

4)  CW  about  E^,  0^,  so  that  Z^  is  vertical  and  is  directed 
north  at  the  radar  site. 

The  transformation  equations  to  accomplish  these  rotations  are 
as  follows: 


‘E4 

*1 

0 

0 

C08X 

s 

0 

-sinX 

s 

- 

0 

COS0 

s 

-sin0s 

0 

1 

0 

.Z4 

0 

- 

sin0 

s 

COS0 

8_ 

sinX 

s 

0 

C08X 

s 

X 


cosX 

o 

0 

- 

sinX 

0 

‘  1  0 

0 

E 

0 

1 

0 

0  cos0 
o 

sin0 

o 

N 

-sinX 

o 

• 

0 

cosX 

0 

0  -sin0 
o 

b 

CO80 

o 

z 

■»  m 
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A)  ORIGINAL  (E.N.Z)  COORDINATES 


C)  2nd  ROTATION:  CCW  ABOUT  N,,X0 
2,  LIES  ON  PRIME  MERIDIAN 


E)  4fh  ROTATION: 
Z,  IS  VERTICAL 


B)  l»*  ROTATION:  CCW  ABOUT  E.  4>0 

N,  ALIGNED  WITH  NORTH  AXIS  OF  E  ARTH 


01  3rd  ROTATION  CW  ABOUT  s 
Z j  LIES  AT  LONGITUDE  X, 


ABOUT 

D  n4  is  directed  north  at  the  radar 


ROTATIONS  TO  RADAR  SITE  COORDINATES 
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To  simplify  the  process  Into  a  product  of  2  matrices  Instead  of 

4,  let  A(0  , A  )  equal  the  product  of  the  first  two  matrices: 
s  s 


AUS,V  • 


1 

0 

- 

0 

cosA 

s 

0 

-sinA 

s 

0 

COS0 

s 

-sin0 

s 

0 

1 

0 

0 

sin0 

s 

cos0 

s 

sin  A 

s 

0 

cosA 

s„ 

A(«  ,X  ) 
s  s 


cosA 


-sin A  sin0  cost) 
s  s  s 

cost)  sin  A  sin0 
s  s  s 


-sin  A 


-sin0  cos  A 
s  s 

cos0  cos  A 

s  9 


Let  B(0  , A  )  equal  the  product  of  the  second  two  matrices: 
o  o 

1 


cos  A  -sin  A  sin0  cos0  sin  A 

o  o  o  o  o 


COS0 


sinf) 


-sinA  -sin0  cosA  cos0  cosA 
o  o  o  o  o 


Notice  that  if  we  ignore  the  subscripts  B(0,A)  “  Aa(0,A). 
The  equation  can  then  be  written: 


r  i 
E4 

N4 

. 

A(0  .A  >  AT (0  » A  ) 
s  s  O  0 

- 1 

Z  PI 

• 

Z4 

TT1 

1 

To  transform  the  vector  CB^  into  the  (E^ ,  N^t  Z^)  coordinate  system: 

L 


CB' 


N 


a(08,a8)  a\0o,ao) 


4U/(4  +  UZ  +  V?) 
4V/(4+u;+V^) 

(4-lT  -  0/(4  +  IT  +  ) 


19 


The  azimuth  of  the  aircraft  is  calculated  hv: 


6  -  ARCTAN 

To  calculate  slant  range,  the  coordinate  system  must  he  re¬ 
scaled,  referenced  to  a  sphere  whose  radius  is  that  of  the  earth. 
A  vector  from  the  center  of  the  sphere  to  the  radar  is  described 
as  follows: 


R 


<E 


0 

0 

+ 


V 


where: 

E  »  Earth  radius  at  the  radar 
s 

Hr  ■  Height  of  the  radar  antenna 
A  vector  to  the  aircraft  is  found  from  CB'^  by: 

—  -  (Es*HT)-ceN  -  (E,  *  V 


where: 

Ht  ■  Aircraft  Height 

Slant  range  p  is  the  difference  between  the  R  and  T  vectors 
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3.0  CONCLUSION 

This  report  explains  a  formulation  used  for  simulated  radar 

data  generation.  The  algorithm  would  be  absolutely  accurate  if  the 

earth  was  a  sphere.  Slight  errors  may  be  introduced  in  conversion 

from  an  ellipsoidal  to  a  spherical  earth  due  to  approximations  in  the 

conformal  latitude  calculation  and  the  calculation  of  earth  radius 

(E  ).  To  minimize  these  errors,  the  conformal  latitude  and  earth 
s 

radius  approximations  should  be  as  accurate  as  possible. 
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